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We present a parameter-free model that estimates the density of excitons in CU2O, related to experi- 
ments that have tried to create an excitonic Bose-Einstein condensate. Our study demonstrates that 
the triplet-state excitons move along adiabats and obey classical statistics, while the singlet-state 
excitons are a possible candidate for forming a Bose-Einstein condensate. Finally we show that the 
results of this study do not change qualitatively in a two-dimensional exciton gas, which can be 
realized in a quantum well. 



PACS numbers: 05.30.Jp, 03.75.Fi, 71.35.Lk, 71.35.-y 



I. INTRODUCTION 

The field of Bose-Einstein condensation jjj has at- 
tracted much attention in recent years, specially in the 
context of trapped alkali- metal atoms j2|. Much ex- 
perimental effort has also been put in creating a Bose- 
Einstein condensate 10] of excitons both in three- 
dimensional crystals , and in quantum wells Jl0| . 

Excitons are bound states of electrons and holes, which 
form in semiconductors after the electrons get excited 
from the valence to the conduction band with use of some 
external probe, a laser beam of light, for example. As 
long as the mean exciton-exciton spacing, which is of 
order n -1 / 3 , with n being the exciton density, is much 
larger than the exciton Bohr radius ag, one can treat the 
excitons as point-like particles that behave as bosons. On 
the other hand for the exciton gas to form a Bose-Einstein 
condensate, the density has to be sufficiently high and the 
exciton temperature T has to be sufficiently low, in order 
for the thermal de Broglie wavelength of the excitons 
to be comparable to the mean exciton-exciton spacing. 
This condition gives rise to the formula connecting the 
critical density n c for Bose-Einstein condensation to the 
critical temperature T c , which in three dimensions and 
for a homogeneous ideal gas gives, 



n c = g{(3/2) 



rnks 
2TTH 2 



3/2 



7^3/2 
1 c 1 



(1) 



where g is the (spin) degeneracy, and m is the exciton 
mass. 

CU2O is a material often used to perform such exper- 
iments with excitons JH|— |^] . The goal is, for typical tem- 
peratures determined by the lattice temperature (i.e., 
liquid helium temperatures), or even somewhat higher, 
to achieve densities that are in the range we described 
above: high enough for the quantum effects to show up, 
but low enough so that the excitons still behave as point- 
like particles, and therefore as bosons. In terms of real 
numbers, assuming, for example, that the temperature is 
T = 2 K, the corresponding density has to be of order 
10 cm~ 3 . Under such conditions the mean exciton- 
exciton spacing is a few hundreds angstroms, which is 



much larger than the exciton Bohr radius (as ~ 5.29 A 
in Cu 2 0). 

In these experiments spectroscopic methods have im- 
plied that the spin-triplet orthoexcitons are highly degen- 
erate [||-(7) moving parallel to the phase boundary. More 
specifically their chemical potential which is extracted by 
fitting the phonon-assisted recombination lines to Bose- 
Einstein distributions is found to be very close to zero. In 
addition the same method implies that the spin-singlet 
paraexcitons have formed a Bose-Einstein condensate [0. 
The lines parallel to the phase boundary given by Eq. (u\) 
are adiabats, i.e., along them the entropy per exciton 
is constant, which is crucial in the discussion that fol- 
lows below. However, more recent studies || of the same 
system have estimated the number of excitons created 
and the volume they occupy, and have shown a clear dis- 
crepancy, since the densities extracted are two to three 
orders of magnitude lower; for such low densities, even 
at the lowest temperatures, the exciton statistics should 
be completely classical. 

In the present study we give theoretical estimates for 
the density of the ortho and paraexcitons as function 
of their temperature, based on a parameter-free model. 
We also give suggestions and predictions for future ex- 
perimental studies, which would put an end to the dis- 
crepancies described above. We also examine the con- 
ditions under which a Bose-Einstein condensate would 
form, since, as we show, this is possible for the paraex- 
citons. Finally, in connection with excitons in quantum 
wells, we demonstrate that our conclusions are unaffected 
for a two-dimensional exciton gas. 

This paper is organized in the following way: we give 
in Sec. II a simple and analytical argument for the fact 
that the orthoexcitons move along adiabats We 
estimate the density of the orthoexcitons and from that 
we extract the paraexciton density in Sec. Ill, discussing 
the long-time behavior of the paraexcitons. In Sec. IV we 
analyze the absorption spectrum of infrared electromag- 
netic radiation inducing transitions from the Is to the 
2p level JlJ] under realistic conditions, and in Sec. V we 
investigate the behavior of a quasi-two-dimensional exci- 
ton gas. Finally we examine the two dominant mecha- 
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nisms that determine the behavior of the orthoexcitons 
in Sec. VI, and summarize our results in Sec. VII. 
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II. ADIABATIC BEHAVIOR OF THE 
ORTHOEXCITONS 

Let us first of all assume that the orthoexcitons and 
the paraexcitons occupy the same volume inside the crys- 
tal, and that they are described by a common tempera- 
ture. Our claim is that the orthoexcitons in this system 
move along adiabats as a result of the balance between 
two competing mechanisms: a heating mechanism that 
comes from an ortho-plus-ortho to para-plus-para con- 
version mechanism fll3|] , and an acoustic-phonon cooling 
process |Q. A similar argument was used in Ref. p4[ , 
but the heating process was assumed to be an Auger 
mechanism ]l5| , |l6| ]. However, in view of the more recent 
experimental measurements || and theoretical calcula- 
tions this Auger mechanism has a negligible rate. 
The model we present here is essentially free of adjustable 
parameters, which makes it very reliable 17 1. 

Let us therefore give the general argument, and then 
examine each mechanism separately. We know from ele- 
mentary thermodynamics that 



dU = TdS - p dV + fi dN , 



(2) 



where Ui is the internal energy, Si is the entropy, pi is 
the pressure, fii is chemical potential, Ni is the number 
of species i, and V is the volume of the exciton gas. 

Solving Eq. (||) in terms of TdS , and expanding the 
differential d(S /N ) we get 
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We show below that the orthoexcitons are in the classical 
regime, and for an ideal monatomic classical gas \i a + 
TS /N = 5k B T/2, and also Po V = N Q k B T. Therefore, 
Eq. (|^) can be written as 
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The differential dll is determined by the two mechanisms 
we mentioned. More specifically, as discussed in Sec. VI, 
the heating rate is equal to an Q , a being a constant, i.e., it 
scales linearly with the orthoexciton density n a , whereas 
the acoustic-phonon cooling rate is equal to —bT 3 / 2 (for 
T 3> Ti, where T) is the lattice temperature), b being a 
constant, i.e., it scales as T 3 ! 2 . Therefore Eq. (Eh implies 
that 
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Let us now examine the terms on the right side of Eq. (|^). 
It turns out that the first two are the dominant ones. 
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FIG. 1. The orthoexciton trajectories on the den- 
sity-temperature plane. The solid line gives the root of the 
polynomial on the right side of Eq. (||), whereas the dashed 
line gives the root of the polynomial on the right side of 
Eq. (f 



More specifically as we show in Sec. VI the order of 
magnitude of these two terms (i.e., the typical energy- 
exchange rates that come from the phonon cooling and 
the ortho-to-para exchange heating) is > 10 meV/ns. 
On the other hand, concerning the last term that comes 
from the expansion of the orthoexcitons, a detailed study 
of this problem has given that the diffusion constant is 
D ~ 10 3 cm 2 /s pOfl . Even at the initial stage of the 
expansion of the cloud when its radius is R > 10 fim, 
dV/(Vdt) - D/R 2 < 1 ns- 1 ; thus this term is < 1 
meV/ns for T « 10 K. Concerning the third term on the 
rig ht side of Eq. (§) one has to distinguish between two 
different cases corresponding to two different experimen- 
tal conditions: For long-pulse laser excitation, of order 
of a few tens of nanoseconds, the population of orthoex- 
citons is essentially determined by the laser and it fol- 
lows its profile (see, for example Fig. 4(a) of Ref. ||). 
Therefore the characteristic timescale tl that is related 
to this mechanism is a few tens of nanoseconds, and the 
corresponding term entering Eq. (EJ) is ~ k B T/rL, or < 
0.1 meV/ns, and thus negligible. Under these conditions 
Eq. (||) takes the form 
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Equation (g) is one of the basic results of this study. To 
see what this equation implies, it is instructive to give the 
two following formulas connecting the entropy, number, 
density, and temperature of the orthoexcitons (which are 
valid for any ideal Bose gas) jH) 
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with r(n) being the gamma function. From the above 
equations one sees that along lines of constant [i /k B T, 
i.e., of constant z , S /N and n /T 3 / 2 are both con- 
stant. Therefore along adiabats, i.e., along lines with 
constant S /N , n a cx T 3 / 2 and according to Eq. (|^), the 
orthoexcitons approach an adiabat, and more specifically 
the one for which n Q = {a/b)T 3 / 2 . This is a line of "equi- 
librium" for the orthoexcitons, which comes as a balance 
between the heating and the cooling mechanisms. Using 
the values for the constants a and b given in Sec. VI, this 
adiabat is [see the solid line in Fig. 1] 



n o (10 
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where the notation n o (10 16 cm -3 ) means that the density 
is to be measured in units of 10 16 cm -3 , and correspond- 
ingly the temperature in Kelvin. The coefficient 10~ 2 in 
Eq. ( |l"o| ) should be compared with w 3.3 given by Eq. (Q), 
which gives the phase boundary, and therefore this adia- 
bat is in the regime where the excitons should not exhibit 
any sign of quantum degeneracy, since —fi /k B T ^ 1. 
For a typical exciton temperature of order 20 K, the 
above equation implies that the orthoexciton density is of 
order 10 16 cm -3 , which is consistent with that estimated 
experimentally in Refs. 0. 

One can also study the dynamics of orthoexcitons. To 
do this, it is convenient to introduce the dimensionless 
quantity a = — ^t /fcsT and rewrite Eq. (13) in the form 
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where a* is the value of a along the specific adiabat with 
n = {a/b)T 3 / 2 . Also 
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Equation ([LlJ) has to be solved along with the rate equa- 
tions describing the change in the exciton density of each 
species, 
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where Gi(t) is the laser production rate of excitons, Tu 
are the intrinsic radiative lifetimes, is the rate for the 
ortho-to-para conversion process, and t~1 is the rate for 
the (reverse) para-to-ortho conversion process. As shown 
in Ref. |jl3), t~* rj cn Q (c is a constant that is given 
in Sec. VI), whereas as explained in Sec. VI the reverse 
process is thermally supressed, r~„ <C t~^. In addition 
the radiative lifetimes of excitons are very long (longer 
than microseconds), and thus Eqs. ( |l4| ) and (15) take the 
simple form 



dflo ^ q /^x tl q 

dt To,p 
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~dt 
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Equations ( p^| ) and (p^|), viewed as densities being func- 
tions of time, do not depend on the exciton temperature 
and can be integrated. Therefore Eqs. (p"l|), (plf), and ( |l6| ) 
can be integrated to give n (t) 7 n p (t), and T(t) [0. How- 
ever we are interested in n G (T), i.e., in the dependence 
of n D on the exciton temperature T, and that is given di- 
rectly by Eq. (|ll|), which implies that the orthoexcitons 
approach the adiabat with a — a» on a timescale of order 
r„ , which is of order of a few tens of picoseconds under 
typical exciton temperatures. 

In another class of experiments short laser pulses (on 
the order of 100 ps) have been used in order to excite 
the crystal. In this case one must include the third 
term on the right side of Eq. (|^) which, however, does 
not change appreciably our basic results. More precisely 
Eq. (||) takes the following form in this case 



dt \ N n 



hT m _ ^k B Tn D . (17) 



The dashed line in Fig. 1 shows the root of the polynomial 
an a - bT 3 / 2 - {hc/2)k B Tn = 0. Although the deviation 
from the adiabatic behavior is not substantial, it is in 
agreement with the experimental observations, as one can 
see in Fig. 2 of Ref. [Q, for example. The open circles 
in this graph correspond to short laser pulse excitation 
(~ 100 ps) and they show some small deviation from the 
adiabats, as opposed to the solid circles (corresponding to 
long pulses, ~ 10 ns) that follow very closely an adiabat. 

From the discussion we have presented up to now there 
is a contradiction: on the one hand spectroscopically the 
orthoexcitons show a high degree of quantum degeneracy 
On the other hand, we argued that the orthoex- 
citons should not show any sign of quantum degeneracy 
(under equilibrium conditions), in agreement with the 
experimental data presented in Refs. |)|. We will not 
address this issue here, which is an open question for 
future studies. 
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III. QUANTUM DEGENERACY OF 
PARAEXCITONS 

Let us now turn to the paraexcitons. Their radiative 
lifetime has been determined to be on the order of mil- 
liseconds ||,|2^], which allows them to form a cold gas, 
with a temperature very close to XJ on these timescales. 
Since on such timescales and for such low temperatures, 
all the excitons will have converted to paraexcitons, they 
are also expected to establish a relatively high density, 
which we estimate now, given the orthoexciton density 
we got earlier. Statistically due to the multiplicity of the 
orthoexcitons, one expects that N a /N p = 3. Therefore 
the paraexciton density (assuming that all the orthoex- 
citons get converted into paraexcitons) is four times the 
one we found earlier. However, as we argued the temper- 
ature is expected to get very close to that of the lattice 
for late times, i.e., for times much larger than the exci- 
ton - acoustic phonon scattering time (on these timescales 
the ortho-to-para conversion process will have converted 
all the orthoexcitons to paraexcitons and there will be 
no heating due to this mechanism). Therefore for an 
initial exciton temperature Tj = 40 K, the paraexciton 
density can get up to 10 17 lf /2 (40 K) cm" 3 , as Eq. @ 
implies, where Ti is to be measured in units of 40 K. 
Even for a density of 10 17 cm -3 the critical temperature 
for Bose-Einstein condensation of the paraexcitons is ap- 
proximately 2 K, i.e., liquid- helium temperatures. One 
important conclusion of this analysis is, therefore, that 
the paraexcitons in the experiments that have been per- 
formed should be very close to the phase boundary, 
or they might have even crossed it already. In the ex- 
periment of Ref. Q the paraexcitons were reported to 
have crossed the phase boundary for Bose-Einstein con- 
densation. The whole analysis was based on spectroscopi- 
cally analyzing the (only) phonon-assisted recombination 
line of paraexcitons, which is very weak, and is close to 
other, much stronger lines. An alternative way of prob- 
ing experimentally the degree of quantum degeneracy of 
paraexcitons is presented in the following section. 

IV. ABSORPTION SPECTRUM OF RADIATION 
INDUCING INTERNAL TRANSITIONS: 
PROBING THE QUANTUM DEGENERACY 

In a possible experiment |p2| that has been proposed 
and has been studied theoretically in Ref. |l2|, it was 
shown that the absorption spectrum of infrared radia- 
tion inducing internal transitions of the excitons from 
the Is to the 2p level is very sensitive to the degree of 
quantum degeneracy of the gas. Therefore, in such an ex- 
periment one should be able to observe the contribution 
of the orthoexcitons and the paraexcitons to the absorp- 
tion separately, with an energy separation of order AE 
(assuming that the width of each distribution is of order 
ksT <C AE). Here AE is the energy splitting between 



the orthoexcitons (lying higher than) the paraexcitons 
due to the exchange interaction fyl. At the zone center 
AE « 12 meV in CU2O, which corresponds to approx- 
imately 150 K. In a temporal study of this experiment, 
following a short laser-pulse excitation (like the one al- 
ready used of a few hundred picoseconds) the contribu- 
tion of the orthoexcitons to the absorption of infrared 
radiation would vanish within a timescale of order T o p , 
i.e., of order nanoseconds under typical conditions. On 
the other hand the contribution of the paraexcitons to the 
absorption would last much longer, on a timescale of or- 
der of their radiative lifetime, i.e., milliseconds p^2|]. As 
shown in Ref. jl2| the appearance of two distinct peaks 
in the absorption spectrum of paraexcitons would signal 
the presence of a Bose-Einstein condensate, since in the 
condensed phase one deals with a two-component system 
and the two peaks would indicate the two different col- 
lective modes of it. But even if the paraexcitons have not 
crossed the phase boundary, but they are highly degen- 
erate, that would still show up clearly in the absorption 
spectrum |L2| . An additional advantage of this method is 
that it does not depend on the strength of the radiative 
recombination lines (which is very weak for the paraex- 
citons, as we mentioned in the preceding section). 

It is interesting that if the lattice temperature increases 
and becomes of order the orthoexciton-paraexciton split- 
ting, that would result in more or less equal rates for the 
ortho-to-para and the para-to-ortho conversion processes, 
and would decrease the paraexciton density substantially 
|p3| . In such a case the orthoexciton lifetime would be 
determined by the radiative lifetime, i.e., microseconds 
P,p2|. Therefore keeping T/ as low as possible is very 
crucial, since it (i) enhances the paraexciton density, and 
(ii) determines the temperature of the paraexcitons for 
late times. 



V. A TWO-DIMENSIONAL EXCITON GAS 

As a final remark we comment on the possibility of con- 
fining the excitons and creating a quasi-two-dimensional 
gas, like in quantum wells |ic| |. As we show in Sec. VI in 
two dimensions the acoustic-phonon cooling rate scales 
linearly with the temperature for T 3> XJ. On the other 
hand the ortho-to-para conversion mechanism is unaf- 
fected by the dimensionality of the system. Thus the 
argument we gave earlier in Eq. (^) now takes the form 

T 7 t {w)»°°°- «■• < 18 > 

where a is the orthoexciton surface density. For S /N 
to be constant in two dimensions a a oc T, and therefore 
according to Eq. ( |l8| ) the orthoexcitons in CU2O would 
still move along an adiabat a a = (b/a)T. In addition a 
two-dimensional gas is expected to undergo a Kosterlitz- 
Thouless transition |24|] to a superfluid state along lines 
on which a a tx T and thus the orthoexcitons would 
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also move along the phase boundary for the Kosterlitz- 
Thouless transition without crossing it, as opposed to the 
paraexcitons that would have the chance of undergoing 
this phase transition. 

VI. ORTHO TO PARA CONVERSION AND 
ACOUSTIC-PHONON COOLING 

Let us now examine in detail the two processes we 
mentioned earlier, starting with the spin-exchange mech- 
anism, where two orthoexcitons collide, exchange their 
electrons or their holes, resulting into two paraexcitons 
in the final states. This process has been studied in detail 
in Ref. |u| and the decay time r 0iP has been calculated 
to be in very good agreement with experiment with 
essentially no adjustable parameter. The decay rate was 
shown to be given by 

t~1 k, 5n o (10 16 cm~ 3 ) ns _1 . (19) 

The reverse process is thermally supressed for tempera- 
tures much lower than AE/ks- The ortho-to-para con- 
version is also a heating mechanism, since there is a gain 
of energy AE per particle. Therefore the heating rate 
due to this process is given by 

1 fdUA = M s , 60 n o (10 16 cm- 3 )meV/ns. (20) 

N o \ dt J o,p T o,P 

We conclude with the cooling process of excitons due to 
their scattering with phonons. The excitons interact with 
the lattice vibrations, and since the lattice temperature 
Ti is kept low (i.e., a few K), the excitons, which in gen- 
eral have a higher temperature, are cooled down. At the 
low energies we consider here the process is dominated 
by collisions with acoustic phonons. The corresponding 
cooling rate per particle is given by Q 

^(^) P ^- a56r3,2(K) ( 1 -?) meV/ " 8 - <21) 

The above cooling rate scales as T 3 / 2 for T ^> Tj, with 
the matrix element due to deformation potential theory, 
the acoustic-phonon energy, and the density of states in 
three dimensions each contributing a factor of T 1 / 2 (they 
scale linearly with the momentum exchange) . For a two- 
dimensional gas the density of states is constant, and the 
cooling rate is proportional to T is this case. 

VII. SUMMARY 

In conclusion, based on a parameter-free model, we 
have given estimates for the density of the excitons in 
CU2O as a function of temperature. The triplet state ex- 
citons are expected to move along adiabats, and to obey 
classical statistics. We have argued that this adiabatic 



behavior is a result of the competition between a spin- 
exchange heating mechanism, and an acoustic-phonon 
cooling process. An open question that needs to be inves- 
tigated is the chemical potential of orthoexcitons, which, 
although they obey classical statistics, is very close to 
zero. 

In addition, according to our model, the singlet-state 
excitons, are a possible candidate for forming a Bose- 
Einstein condensate. Finally, we have shown that in a 
two-dimensional exciton gas the basic conclusions of our 
study remain the same. 
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